Coupled electron and phonon transport in one-dimensional atomic junctions 
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Employing the nonequilibrium Green's function method, we develop a fully quantum mechanical 
model to study the coupled electron-phonon transport in one- dimensional atomic junctions in the 
presence of a weak electron-phonon interaction. This model enables us to study the electronic and 
phononic transport on an equal footing. We derive the electrical and energy currents of the coupled 
electron-phonon system and the energy exchange between them. As an application, we study the 
heat dissipation in current carrying atomic junctions within the self-consistent Born approximation, 
which guarantees energy current conservation. We find that the inclusion of phonon transport is 
important in determining the heat dissipation and temperature change of the atomic junctions. 
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I. INTRODUCTION 



The electronic transport and phononic transport in 
meso- and nano-structures have attracted a great deal 
of interest in the past two decades, although their devel- 
opment is not so parallel sometimes. These structures 
display important quantum effects due to the confine- 
ment in one or more directions^. The quantized electrical 
conductance* was observed much earlier than that of the 
thermal conductance^ mainly due to the difficulty in mea- 
suring the thermal transport properties. Electrons and 
phonons are not two isolated systems. Their interactions 
are important for both electronic and phononic trans- 
port. With the development of both fields there arises 
the requirement to study the coupled electron-phonon 
transport from time to time. When studying electronic 
transport problems, one usually assumes that electrons 
interact with some phonon bath where the phonons are in 
their thermal equilibrium state characterized by the Bose 
distribution. This simple assumption is not able to give 
satisfactory results in some cases where the phonons are 
driven out of equilibrium by the electrons. This is espe- 
cially true in places where the thermal conductance is low 
or the phonon relaxation is slow^^. To take into account 
the nonequilibrium phonon effect, one usually introduces 
into the electronic transport formalism some phenomeno- 
logical parameters that describe the phonon relaxation 
process. In engineering applications, as the size of the 
electronic devices decreases to nanoscale, the heat dissi- 
pation and conduction in these structures become crit- 
ical issues, which may influence the electronic proper- 
ties dramatically- 6 -. Only studying the electronic trans- 
port is not enough in these cases. On the other hand, 
heat transport in one-dimensional (ID) structures has re- 
ceived considerable attention recently 6,7,8 . Fourier's law 
of heat conduction is no longer valid in many ID systems. 
The microscopic origins of the macroscopic Fourier's law 
remain one of the most frustrating problems in nonequi- 
librium statistical mechanics. Since the electrons and 
phonons both contribute to the heat conduction, their 
relative roles in many nanostructures are still not clear. 



Especially in semiconductors, which one carries the ma- 
jority of the thermal current is not a trivial problem. To 
answer these questions, we need some general models, 
which take into account the electron, phonon transport, 
and their mutual interactions. 

Theoretically, although the development of electronic 
transport in ID structures has been very striking, that 
of the phononic transport is relatively slow. Classical 
molecular dynamics (MD) and the Boltzmann-Peierls 
equation are the widely used methods in phononic trans- 
port. MD method is not accurate below the Debye 
temperature, while the Boltzmann-Peierls equation can 
not be used in nanostructures without translational in- 
variance. In both cases, the quantum effect becomes 
important 1 . Only recently, the nonequilibrium Green's 
function mcthod^ ' 10 i 1:L i 12 ; which has been widely used 
to study the electronic transport, has been applied to 
study the quantum phononic transpor t 13 i 14 ' 15 i 16 i 17 . As 
far as we know, the study of the coupled electronic and 
phononic transport in nanostructures is rar o 18 ' 19 i 20 i 21 . 
In Ref. [H the authors considered the nonequilibrium 
phonons in molecular transport junctions. Galperin and 
co-authors analyzed the heat generation and conduction 
in molecular systems^. In this paper, using the nonequi- 
librium Green's function method, we study the coupled 
electronic and phononic transport in ID atomic junc- 
tions. The formalism is similar to that of Ref. [HI In our 
model the electron subsystem is described by a single- 
orbital tight-binding Hamiltonian, and the phonon sub- 
system is described in a harmonic approximation. We 
assume that the electron-phonon interaction is weak so 
that the perturbative treatment is valid. The strong- 
interaction case is the scope of future work. 

The rest of the paper is organized as follows. In Sec.HTl 
we introduce the ID model system, and derive expres- 
sions for the electrical, energy current of the coupled 
electron-phonon system. In Sec. IIIII we show the heat 
generation in one- and two-atom structures under differ- 
ent model parameters. Sec. IIVI is the conclusion. In Ap- 
pendix [X][C] we give some technical details of our deriva- 
tion. 
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II. COUPLED ELECTRONIC AND PHONONIC 
TRANSPORT 

A. The Hamitonian 

Our model system is an infinite ID atomic chain as 
shown in Fig. [T] The electrons and atoms are only al- 
lowed to move in the longitudinal direction. We treat 
the atoms as coupled harmonic oscillators, and take into 
account their nearest neighbour interactions up to the 
second order. We assume that there is only one single 
electronic state for each atom and take into account hop- 
ping transitions between the nearest states. This corre- 
sponds to a single-orbital tight-binding model. Also, we 
assume that there is only one spin state for each orbital. 
Following Caroli^, we divide the whole system into one 
central region and two semi-infinite leads, which act as 
electrical and thermal baths (Fig. [TJ. The Hamiltonian 
of the whole system is 

H = E H P 

a=L,C,R;0=e,ph 

+ Y, (Hf + H^)+H cph . (1) 

a=L,R;/3—e,ph 

The electron-phonon interaction Hamiltonian -ff Gp h is 
non-zero only in the central region. The electron Hamil- 
tonian reads 

H« = Y,£?4 a c?+ E *S4 a # ( 2 ) 

i |* — J 1=1 

where cJ Q and cf are the electron creation and annihila- 
tion operators, ef is the electron onsite energy, and tfj is 
the hopping energy between adjacent states, i and j run 
over the sites in the a region. The coupling Hamiltonian 
with the leads is 



H. 



LC 



LC JL C 
ij L i <-j i 



(3) 



and 




FIG. 1: Shematic diagram of the ID coupled electron-phonon 
system and the parameters used in the model. The big dots 
in the bottom line represent atoms, while the small dots in 
the upper line represent electron states. They are coupled via 
the electron-phonon interaction. 



— Kq / y/mfm/j (i ^ j). Here Kq is the spring constant, 
and mf is the mass of the ith atom in the a region. Like 
the electrons, the coupling Hamiltonian with the leads is 



ttLC _ 

H ph - 



and 



ttCR _ 1 

^ P h - - 



ij 



(6) 



(7) 



We also have K Ca = K aC \ 



The non-zero elements are 



and Kgg+v 



K LC K CL K RC 
-"■01 ' -"-10 i -"-n+l,! 

The electron-phonon interaction is included within the 
adiabatic Born-Oppcnhcimcr approximation. First, the 
electron subsystem is solved with all the atoms in their 
equilibrium positions. Then, the isolated phonon sub- 
system is considered. After that, the electron-phonon in- 
teraction is turned on by allowing the atoms to oscillate 
around their equilibrium positions. Within this picture, 
the electron-phonon interaction isii 



cph 



5> 

i,j,k 



k c\cjU k . 



The interaction matrix element is M*j = ( / 



dH,. 



du k 



(8) 



All 



the operators in Eq. ((SJ) are in the central region, so we 
omitted the superscript C. In our model, the electron 
operators are in the second quantization, while that of 
the phonons are in the first quantization. 



H CR = 



(4) 



B. Green's functions 



H^ L and H^ c have similar expressions. We also have 
t aC = t c a ^ a = L R For Qur 1D tight-binding model, 
t has only one non-zero element. If we label the central 
atoms with indices 1 to n as shown in Fig.[Tl the non-zero 
elements will be tj£f, tf L , , 
The phonon Hamiltonian is 



an d *n,n+l- 



Ha 
ph 



^E« 



E 

-3 1=0 



n 3 



(5) 



uf and iif are the mass-renormalized atom displacement 
and momentum operator. K% = 2Kq /mf, and Kf^ = 



The nonequilibrium Green's function method for the 
electronic transport is discussed in Refs. |1 10lllllll2l a nd 
that for the phononic transport in Refs. l^I33lilI3- 
Here we concentrate on the electron-phonon interactions. 
The definition of the electron contour-ordered Green's 
function is Gjfe(r, r') = -i{T{cj(r)cl( , r')}) ) and the 
phonon counterpart is Djk(r, r') = —i(T{uj(T)uk(r')}) . 
Here t is time on the Keldysh contour, and 7~{- • •} is 
the contour-ordered operator. We set ft = 1 through- 
out the formulas. Without the electron-phonon inter- 
action, the isolated electron and phonon problem can 
be solved exactly. We denote these Green's functions 
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as C?o(t,t') and Dq(t,t'), respectively. In our case, 
it is convenient to write the Hamiltonians as matri- 
ces and work in the energy space. The electron re- 
tarded and advanced Green's functions are G r (e) — 
Ggt(e) = [{e + iri)I-H^ -W L {e)-W R {e)]-\ lis 
an identity matrix, and r\ — > + . The retarded self- 
energy Y7 a — t Ca g 7 a t aC is due to the interactions with 
the lead a. The retarded Green's function of the semi- 
infinite lead g r a can be obtained analytically (Appendix 
|A"|) . The "less than" Green's function is given by Gq = 
G5(E< + £<)G§, where £< = -£(E£ - E£). /• is 
the Fermi-Dirac distribution. The phonon retarded and 
advanced Green's functions are2 3 Dq(uj) = Dq'(cj) = 

[(u + iri) 2 I- K c -W L (u})-W R (u)]~ l . The lead re- 
tarded self-energy is U r a {cj) = K Ca d r a (uj)K aC . d r a also 
has analytical expression (Appendix lA"]) . The phonon 
"less than" Green's function is D< = D r (U< + n<)Dg, 
where IT< = fP h (Il r a -II"). /P h is the Bose distribution 
function. 

Knowing the bare electron and phonon Green's func- 
tions Go and D$, we can include their interaction as per- 
turbation. Following the standard procedure of nonequi- 
librium Green's function method, we can express this 
interaction as self-energies. The full Green's functions 
are obtained from the Dyson equation, e.g., for elec- 
trons G r ' a = Go 



y< 



y< 



GJX'phG^, and G< 



G r S<G a . 



is the total self-energy. Keeping 
the lowest non-zero order (the second order) of the self- 
energies, we have two (Hartree- and Fock-like) terms for 
the electrons, and one polarization term for the phonons. 
This is the so-called Born approximation (BA)ii. The 
Fock self-energies are 

G< H (e-cj)D< kl (cj)^-M l in , (9) 



££<( E )= iM* 



and 



1(e) 



iMt 



du> r 

2^ 



2tt" 



G r 0i Ae-u)D< kl {w) 



+G r 0ij (e-u)D r m (u J ) 



M l . 



(10) 



The "less than" Hartree self-energy is zero, and the re- 
tarded one is 

s = -iMLWuV = mil J G^ lk (s)^. (ii) 

This term is a constant for all energies, which represents 
a static potential due to the presence of phonons. The 
self-energies for the phonons are 



and 



(12) 



+GZ kl (e)G a Je-u) 



M'' 



(13) 



In Eqs. (|9l fl3| . sum over internal indices is assumed. The 
self-consistent Born approximation (SCBA) is obtained 
by replacing all the bare Green's functions Go and Do in 
Eqs. tUBD with the full G and D±I. In Appendix [5J 
we show that the SCBA fulfills the electrical and energy 
current conservation, while BA fails. 



C. The electrical and energy current 

The electrical and energy current can be expressed by 
the Green's functions. The electrical current out of the 
lead a i o 11 ' 24 

f dr 

Ja =e J 27 Tr {G > (£)S<(e)-G<( £ )E>( £ )}. (14) 

The electron energy current is 

Jf c = | JeTr{G>( £ )£<( e )-G<( £ )£>(e)}. (15) 

The electron heat current is obtained from Eqs. (fTH - 
IT5|) as J 1 *' — Ja'° — HaJa/e- Ha is the lead chemical 
potential. The derivation of the phonon energy current 
runs parallel with that of the electrons 1 ^ 

J « Ph = Tr{D>Hn<H - D < (uj)U> (uj)}. 

(16) 

For phonons the energy current is the same as the heat 
current. When there is no electron-phonon interaction, 
the electron energy current is conserved throughout the 
structure. So is the phonon energy current. In the pres- 
ence of such an interaction, only the total energy current 
is conserved due to the energy exchange between them. 
The phonons do not carry charges, so in both cases the 
electrical current is conserved. Since we can't get the 
exact self-energies in most cases, we need some approxi- 
mations. Properly defined self-energies should fulfill the 
electrical and energy current conservation 



a 



(17) 



(18) 



where a runs over all the leads. We justify that the SCBA 
fulfills these conservation laws, while the BA fails to con- 
serve the energy current (Appendix B). Provided we sat- 
isfy these conservation laws, we can write the electrical 
and energy current in symmetric forms. The electrical 
current is 



J = ej ^f c ( £ ) [fl{e) 



f e R (e)}. 



(19) 



The transmission coefficient reads 



G r (rL + -r cp h — S c )G a Tft 



-G r TLG a (Tji + -r cp h 



S e ) 



(20) 



4 



where S e is 



5° 



l(/^ + /£)r eph + *E< h 

fC fC 

J L J R 



(21) 



r Q = — E£), a = L,R is the electron level-width 
function. r ep h = ph — S" ph ) is due to the electron- 
phonon interaction. The total energy current is 



J E = 



1*. 
2 



(22) 



The phonon transmission coefficient is 



f Ph = Tr <! - 



D r (A L + iA cph - S^Z^A* 



-D r A L D a (A R + ±A cph + Sv h ) }, 



where 5 ph is 



i(/£ h + /r)Ae ph -ffl 



5 ph 



< 

eph 



fPh 



(23) 



(24) 



A Q = «(n^, — 11^) is the phonon level- width function. 
A cp h = *(ne P h~n° ph ) is due to the electron-phonon inter- 
action. Eqs. (|19H24j) are the generalization of the Caroli 
formula— to include the electron-phonon interaction. 



III. HEAT GENERATION IN CURRENT 
CARRYING ID ATOMIC JUNCTIONS 

As an application of the formalism in Sec. [ill we 
study the heat dissipation in current-carrying ID atomic 
junction a 18 ' 25 i 26 ' 27 i 28 ' 29 i 30 ' 31 . In the presence of potential 
difference between the two leads, there will be electrical 
current flowing between them. When the electrons pass 
the central region, there is energy exchange between the 
electron and phonon systems. The energy dissipated into 
the phonon system makes the atom temperature higher 
than that of the leads if it is not efficiently conducted to 
the leads. If the electron-phonon interaction is weak, the 
energy dissipated into the phonon system is only a small 
fraction of the electron energy current. But this small 
fraction still influences the transport properties of the 
atomic junction and even leads to junction breakup 32 ' 33 , 
especially when the thermal conductance is low. Dif- 
ferent models have been used to study the local heating 
effect. Some simply assume that the phonons are in their 
thermal equilibrium states 3 -^. Some take into account the 
phonon transport by using the rate equations 3 ^ or other 
semi-classical model o 26 ' 27 ' 32 . Few of them take into ac- 
count the quantum effect in heat transpor t 18 ' 34 . Our 
model treats the electron and phonon transport on an 



equal quantum-mechanical footing, and includes their in- 
teractions self-consistently. The heat generation is given 
by (Eq. JB9])) 



w)M\ n ]. (25) 



[G> m {e)M^D<{uj)G<Ae 



At zero temperature, we can get an analytical expression 
Eq. (jCip for a single-atom structure by using the bare 
Green's functions Go and Dq in Eq. (|25p (Appendix [Cj) . 
Equation (|C1|) can reproduce most qualitative features of 
heat generation in a single atom, except that it does not 
take into account heat conduction in the phonon system. 

We first study the case where the lead energy band is 
wide compared to the voltage applied to the structure. 
For most metallic leads, this condition should hold. In 
the weak electron-phonon coupling regime, the Born ap- 
proximation should give acceptable results for the heat 
generation, although physically it is not a good approx- 
imation. Figure [2] shows the heat generation of a single 
atom (n = 1 in Fig. []} computed using Eqs. (|B5[) and 
(|B8|) under BA and SCBA, respectively. The parameters 
used in the calculation are stated in the figure caption. 
With these parameters, the electron energy band is in 
the range — 1 < e < 1 eV. The chemical potential of each 
lead is zero in equilibrium. The phonon energy is ap- 
proximately lo — 0.05 eV. In all the results presented in 
this section, the temperature is T = 4.2 K, the electron- 
phonon coupling matrix M = 0.08 eV/(A-amu2). The 
cut-off energy of the electron system is 2.1 eV, and the 
phonon system is 0.2 eV. The energy spacing is dis- 
cretized into grids of 1 meV. Equation (|B5[) gives the 
energy decrease of the electron system, while Eq. (|B8[) 
gives the energy increase of the phonon system. Numeri- 
cal results from Eq. (|B5|) and Eq. (|B8j) under SCBA have 
some slight discrepancy. This is due to numerical inac- 
curacies. But most of the discrepancy under BA comes 
from the difference between the bare and the full Green's 
functions, which may become even larger for some pa- 
rameters. So BA should be used with care in the study 
the energy exchange between the electron and phonon 
system. We also note that although Eqs. (|B5[) and (|B7|) 
are equivalent, numerical result from Eq. (|B5[) is unsta- 
ble in many cases. The reason is that the energy ex- 
change between the electron and the phonon system is 
only a small fraction of the total electrical energy cur- 
rent. Equation (|B5P is the difference between two large 
numbers, so our numerical integration has to be accurate 
enough to get a reasonable result— On the contrary, 
Eq. (f2"S"]) is much more stable since we have got the dif- 
ference analytically. All the results presented below use 
this equation. 

From Fig. [21 we can see two threshold values in heat 
generation. The first one corresponds to the onset of 
phonon emission. Under low temperatures, the equilib- 
rium phonon occupation is very small, so the phonon 
absorption process seldom takes place. If the applied 
bias is smaller than the phonon energy, electrons don't 
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0.1 0.2 0.3 0.4 0.5 
Voltage (V) 



FIG. 2: Comparison of different methods to compute the heat 
generation in a single atom structure. The four curves cor- 
respond to results from Eqs. ()B5|I and (|B8[) under BA and 
SCBA, respectively. If we label this single atom as index 1, 
its electronic onsite energy is written as =0.1 eV. The 
onsite energy of the leads is e L = e R = eV. The hopping 
energy is th = tf* = 0.5 eV. The non-zero electronic cou- 
pling with the lead is ff £ = t$f = 0.1 eV. The matrix ele- 
ment of the sing le atom is Ku = 0.654 eV/(A 2 - amu). The 
spring constant between the lead atoms is Kh = K R = 0.654 
eV/(A 2 -amu). The non-zero atomic coupling with the leads 
is Kg L = K$f = 0.127 eV/(A 2 -amu). 



have enough energy to emit one phonon. So the heat 
generation is zero. Once the applied voltage is larger 
than the phonon energy, phonon emission turns on. The 
heat generation increases almost linearly with the applied 
bias (inset of Fig. [2]). This is different from the electrical 
current, which increases smoothly in this regime. The 
second threshold value corresponds to the alignment of 
the left lead chemical potential with the electron onsite 
energy eV = 2sq (positive bias /i; > fi r ). The electron 
transmission is nearly unity above the onsite energy. The 
larger the transmission, the larger the current and the 
heat generation provided that the other parameters re- 
main unchanged. These two threshold behaviours may 
become less obvious when the coupling with the leads 
get stronger. As a result of coupling, the discrete elec- 
tron and phonon density of states (DoS) extends to a 
small energy region around their discrete values. The 
continuous phonon DoS leads to the broadening of the 
first threshold behaviour, while the continuous electron 
DoS is responsible for that of the second. It is smoothed 
out when the coupling is large enough (Fig. [3]). Only 
electrons whose energies are within the broadened energy 
spectrum can tunnel across the central atom. The heat 
generation reaches maximum when the electron states in 
one lead are all occupied in this energy range, while those 
in the other are all empty. 

The electron-lead coupling not only leads to the elec- 
tron level broadening, but it also influences the elec- 
tron tunneling time. The larger this coupling, the less 
time electrons spend in the central region. In Fig. [3J we 
show the heat generation and the atom temperature for a 



single-atom structure under different electronic coupling 
strengths. The definition of temperature is ambiguous 
in nanostructures^. Here we use the method proposed 
in Ref. [ill We can only see one threshold behaviour at 
about 0.2 V, which is smoothed out when the coupling is 
larger than 0.2 eV. The temperature and the heat gener- 
ation show similar trends. The saturate voltage of heat 
generation increases with the strengthen of the electron- 
lead coupling. This is due to the coupling induced atomic 
level broadening. The decrease of the heat generation 
and temperature with increasing electron-leads coupling 
can be easily understood. The larger this coupling, the 
less time electrons spend at the central atom. Since the 
electron-phonon interaction takes place there, the heat 
generation decreases. We also show the heat generation 
as a function of electron-lead coupling in the inset of the 
lower panel. The applied voltage is 0.3 V. On one side, 
when the coupling is too small, few electrons can tun- 
neling through the atom. The heat generation is small. 
On the other, when the coupling is very large, the elec- 
tron tunneling process is too quick for the phonons to 
interact with the electrons. The heat generation is also 
small. It has a maximum value at some moderate cou- 
pling strength. This is different from the electrical cur- 
rent, which increases monotonously with the increase of 
coupling strength. 




Voltage (V) 

FIG. 3: Heat generation Q and the atom temperature T under 
different electron coupling strength fjjf = t^x = 0.05, 0.1, 
0.2, and 0.3 eV, respectively. Other parameters are the same 
with Fig. [5] The inset shows the heat generation as a function 
of electron coupling strength at an applied bias V — 0.3 V. 

The atom-lead coupling determines how well the gen- 
erated heat can be conducted into the surrounding leads. 
One of the important reasons why we are interested in 
the heat generation in nanostructures is that it may leads 
to temperature increase and even structure breakup. To 
study the temperature change, we need to take into ac- 
count not only the heat generation, but also the heat con- 
duction into the leads. In the simplest one-atom struc- 
ture, the heat conductance is mainly determined by the 
atom-lead coupling. Our model includes this intrinsically. 
Figure [3] shows the heat generation and the atom temper- 
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ature as a function of atom-lead coupling. For the heat 
generation, the BA and SCBA results show large differ- 
ence around the resonant position, which corresponds to 
a perfect atomic junction. For the atom temperature, BA 
and SCBA give almost the same results. In the case of 
a perfect junction, the heat generation reaches its max- 
imum value, while the atom temperature is the lowest. 
The reason is that the perfect junction has the best heat 
conductance. When the atom-lead coupling is weak, the 
heat generation is small. But the poor heat conductance 
can still result in a much higher temperature than the 
surrounding leads. We also show the heat conductance 
as a function of atom-lead coupling in the inset of the 
upper panel, which shows a sharp peak at resonance. 




0.2 0.4 0.6 0.8 
K (eV/A 2 amu) 



FIG. 4: Heat generation Q and the atom temperature T as 
a function of the atom-lead coupling Ki L = = K. 

Dashed and dotted lines correspond to results of SCBA and 
BA, respectively. Other parameters are the same with Fig. [2] 
The inset shows the thermal conductance Ki £LS Et function of 
K. The unit is 1 x 10~ 12 W/K. 

In Fig. we show the heat generation as a function 
of electron onsite energy at different biases. We assume 
that we can tune the the onsite energy via a gate voltage. 
When the applied bias is less than the phonon energy, 
there will be no heat generation. When the bias energy is 
slightly larger than the phonon energy and less than 2w, 
there are two energy positions where the heat generation 
is the largest. These two peaks are approximately at 
— Q.beV + u and 0.5eV — uj. They merge into a single one 
at a bias of eV = 2u> until it reaches saturation. After 
that, this peak broadens, and becomes ladders. All these 
behaviour can also be explained by the analytical result 
ofEq. (ED- 

In Fig. [5] we show the heat generation of a two-atom 
structure (n = 2 in Fig. [I}. The central region has two 
identical atoms. Interaction between them leads to two 
discrete energy levels. One is at eV, and the other at 
0.4 eV. When the electrical coupling between the leads 
and the central region is small (0.1 eV), additional to the 
threshold behaviour at eV = u>, there are two ladders 
corresponding to the phonon assisted resonant tunneling 
across the two electrical levels. If the electrical coupling 
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FIG. 5: Heat generation Q as a function of electrical onsite 
energy ef under different biases. From the inner to the outer 
side, the applied biases are V = 0.05, 0.10, 0.15, 0.20, 0.25, 
0.30, 0.35, 0.40, 0.45, and 0.50 V, respectively. 

gets larger (0.2 eV), the two ladders broaden out. Again 
this is attributed to the coupling induced level broad- 
ening. The heat generation for the two-atom structure 
is much larger than that of a single-atom structure. The 
more the electrical levels, the larger the electrical current 
and heat generation. It is worth noting that for multi- 
atom structures the distribution of the electrostatic po- 
tential may influence the results significantly 3 -. In the 
above calculation, we assume that the two electrical lev- 
els don't change with the applied bias, and that we can 
tune their positions via a gate voltage. 
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Voltage (V) 

FIG. 6: Heat generation Q as a function of applied voltage for 
a two-atom structure. The two-atom onsite energy is e c = 0.2 
eV, the hopping energy is t c — 0.1 eV, and the spring con- 
stant is K c = 0.654 eV/(A 2 -amu). The two leads are iden- 
tical. The electron onsite energy is e L = e R — 0, and the 
hopping energy is t L =t R = 0.5 eV. Their spring constants 
are the same as the central region. The non-zero coupling cou- 
plings with the leads are Xf L = K§f = 0.327 eV/(A 2 -amu), 
and t°a = t^f = 0.1 (solid), 0.2 eV (dashed), respectively. 

If one of the metallic leads is replaced by a semicon- 
ductor, there will be some new features in the electrical 
current and the heat generation. In our simple model, 
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we can alternate the electron onsite energies between two 
values to mimic a simple semiconductor (Appendix [A"]) . 
In Fig. [7] we show the heat generation and the electri- 
cal current for such kind of structure. The alternating 
onsite energies of the left lead are —0.1 and —0.2 eV, re- 
spectively. This produces an energy band-gap of 0.1 eV. 
Other parameters are given in the figure caption. We can 
see that there appears negative differential conductivity 
in the current-voltage characteristics due to the semi- 
conductor band-gap. This qualitatively agrees with the 
experimental 3 ^ and first-principle 3 - 7 " studies. The heat 
generation curve is slightly different. Additional to its 
threshold behaviour, the peak and valley positions are 
also different. The electrical current has a peak when 
the chemical potential of the lead is aligned with the 
central electrical level, while the peak of the heat gen- 
eration shifts to the right by one phonon energy. This 
corresponds to the phonon-assisted resonant tunneling. 
The current and heat generation decrease when the sin- 
gle electrical level is within the band-gap of the left lead. 
The peak-to-valley ratio depends on the coupling with 
the semiconductor lead. In the limit of small band-gap 
and large coupling, we recover the metallic lead results. 



ied the heat generation in one- and two-atom structures 
coupling with different leads under a broad range of pa- 
rameters. Especially, we studied the influence of the ther- 
mal transport properties on the heat generation and atom 
temperature of the central region. The results on semi- 
conductor leads agree qualitatively with the experimental 
and first-principle studies. This model can be easily ex- 
tended to study more realistic structures such as molec- 
ular transport junctions and metallic nanowires. The 
electron, phonon Hamiltonian, their interaction and lead- 
coupling matrices can all be obtained from first-principle 
calculation s 17 ! 38 ' 39 . The surface Green's functions for 
bulk leads can be computed by recursive metho d 17 ' 38 . It 
is also possible to include the electron-electron and the 
phonon-phonon interaction s 14 ' 17 . 
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APPENDIX A: 



0.4 0.6 
Voltage (V) 

FIG. 7: Heat generation Q and electrical current J as a 
function of applied voltage for a single-atom structure. The 
left lead is a semiconductor. Its alternating onsite energies 
are —0.2 and —0.1 eV, respectively. The chemical potential 
is fi L = 0.05 eV higher than the conduction band bottom, 
which corresponds to n-type doping. Other parameters are 
K? L = K£ c = 0.4 eV/(A 2 -amu), t? L = tW = 0.1 eV, 
ef = 0.1 eV, Kg = 0.654 eV/(A 2 -amu), e R = -0.05 eV, 
t L = t R = 0.5 eV, and K L = K R = 0.654 eV/(A 2 -amu). 



IV. CONCLUSION 

We studied the coupled electron and phonon transport 
in ID atomic junctions in the weak electron-phonon inter- 
action regime. Base on the nonequilibrium Green's func- 
tion method, we derived the electrical, energy current of 
the coupled electron-phonon system, and the energy ex- 
change between them. We showed that the SCBA con- 
serves the energy current. Using this formalism, we stud- 



SURFACE GREEN'S FUNCTIONS 
OF THE ID LEAD 



In this Appendix, we show that for the ID tight- 
binding model the lead self-energies can be expressed 
analytically^ 7 -. The electron and phonon self-energies are 
similar in their form. Here we take electrons as an exam- 
ple, and give the phonon results directly. We assume that 
the onsite energies of the electrons alternate between e\ 
and £2. The hopping energy is tfj = to- If £1 = £2, we get 
a continuum band. This corresponds to a metallic lead. 
If they are not equal, we get two bands with a band gap. 
We can take the lower as the valence band (VB) , and the 
upper as the conduction band (CB). We use this method 
to mimic a semiconductor lead. In this case, the semi- 
infinite lead has two electron states in each period. In the 
tight-binding model, only the left- (right-) most state of 
the central region is coupled to the left (right) lead. So 
we only need to know the surface Green's function, e.g., 
for the left lead it is go — go - We assume the retarded 
Green's function is 



state 1, 
state 2. 



(Al) 



Putting it into the definition of the retarded Green's func- 
tions [(£ + irfjl — H]g r = I, we have 

- t Q a + (e + ii] - £ 2 )c 2 - toCiA = 0, (A2) 

- hc2 + (e + ir) - £i)ciA - t Q c 2 X = 0. (A3) 
From Eqs. (|A2IIA3[) , we get an equation for A 
(s + irj — £i)(e + ir/- £ 2 )~ 



A^ + 



2- 



t 2 




A + 1 = 0. (A4) 
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The condition that Eq. (|A4|) has travelling wave solutions 
gives the dispersion relation 



(ei+e 2 )-V(ei-e2) 2 + l 



£ 2 < e < 



(ei+e 2 ) + V(£i-£2) 2 + 



< £ < £1 



Ki/r 



(VB), 
(CB). 



(A5) 



We assume ei < e 2 without loss of generality. The energy 
band-gap is e% — e\. If they are equal, the two bands 
merge into one, which corresponds to a metallic lead. 

For the surface Green's function of the left lead, we 
also have 



(e + it] - £i)ci - t c 2 = 1. 
From Eqs. (|A2IA6[) , we get 



50 



ww ( VB )> 



(A6) 



(A7) 



|A| > 1 is one of the roots of Eq. (|A4|) . The surface 
Green's function of the right lead is identical. 

We can also alternate the atom masses to generate a 
phonon band-gap. In our model the mass change will 
modify the renormalized spring constants. The diagonal 
elements of the dynamical matrix will be two alternating 
values K% — 2h\ or 2&2, while the off-diagonal elements 



will be a single value Kf, — —\/k]k 2 , where 



3 = 1- 



If we assume that k 2 > k\, the acoustic band (AB) is 
< uj 2 < 2fci, and the optical band (OB) 2k 2 < uj 2 < 
2(ki + k 2 ). The surface Green's function is 



do - < '■>< (OB), 



(A8) 



(1+A)fcifc 2 



where tt n = (uj 
of 



it]) 2 — 2k n . |A| > 1 is one of the roots 



A 2 



k x k 2 



A + 1 = 0. 



(A9) 



In all the simulation results of present paper, the two 
spring constants are equal (k\ = k 2 ), which correspond 
a single continuum phonon band. The electron onsitc 
energies are also equal {e\ — e 2 ) except in Fig. [3 where 
we set £1 = —0.2 eV and £2 = —0.1 cV to mimic a 
semiconductor lead. 



APPENDIX B: ENERGY CURRENT 
CONSERVATION 

In this Appendix, we justify that the SCBA satisfies 
the energy current conservation. The justification of the 
electrical current conservation is given in the Refsi 40 ' 41 . 
What we need to prove is that 



J. 



E,phN _ 



= 0. 



(Bl) 



The electron part is 

E^E / ^Tr{G>(£)E<(£)-G<(£)E>(£)}. 



Using the important relatio n 39 ' 40 

Tr{G>£<-G<£>} - 0, 

we get 



(B2) 
(B3) 



E' E 



J «' C = "/ £ £ T^^pJe) ~ G<(£)S> h (£)}. 

(B4) 

The Hartree term does not contribute to the current di- 
rectly. It's just like a static potential which only modifies 
the Green's function. Putting the Fock self-energy into 
Eq. (fB4|) . we have 



l 

jn 



-q = E J « ,e 

a 

= -iJ^J ^e[G> m (e)M^D<(uj)G<(e-uj)M 

-G< m (s)M^D> (uj)G> (£ - uj)M l jn ] . (B5) 

Sum over all the indices is assumed. The heat generation 
Q is the energy decrease of the electron system, which 
should also be the energy increase of the phonon system. 
Replacing u) by —uj, using the symmetric properties of 
the phonon Green's functions 1 -^, replacing e by e — uj, 
and finally changing dummy variables, we get 



de 
2^ 
ds 
2^ 



doj 
to' 
doj 
2^ 



£ [G< m (e)M k mi D>{uo)G>(e - uj)M) n 
£ [G< m (e)M^D>(-uj)G>(e + U )M. 
S / ^ E ~ u) [ G n TO (e-w)M^D<( W )G>( e )M. 

(£ - uj) [G> m (e)M^D<(u)G<(s - w)M. 



1 1 

jn} 



ds 
2^ 



duj 
2^ 



1 1 



(B6) 



Putting Eq. (|B6[) back into Eq. (|B5j). we get 

-Q = E J *' e 

a 

= ~ l J^J ^[ G l m (e)M k mt D<^)G<{e-oj)M l jn ] 
* 0. (B7) 
For the phonon energy current we have 

Q = E^ ph 

a 

' !: f ' Lc [ Dl m {u)Mr k GUe)G>{e-uj)Ml } 

(B8) 



2n / 4w 



-d: 



MM^Gl^Gf^-^Ml 
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Following the same procedure as electrons, finally we get 

Q =E J "' ph 



de f dtu 



2tt / 2?r 



u> [G> m (e)M^D<(u;)G<(e - u)M l jn ] 



So we still have 



£ ( J « c + J « ph ) = °- 



(B9) 



(BIO) 



Eqs. (|B7[ IB9|) give the energy exchange between the 
electron and the phonon system, which is also the heat 
generation of the atomic junction. Replacing D < , G < 
by £><, G< in Eq. fBfi, and G>, G< by G>, G< in 
Eq. (|B9[) , we get the results under BA. We can find that 
the energy increase of the phonons does not equal to the 
energy decrease of the electrons under BA. 



APPENDIX C: ANALYTICAL RESULT AT ZERO 
TEMPERATURE 



in their equilibrium states. Under these approximations, 
the heat generation is (assuming eV > luq) 



Q 



-+U1Q 



ds 
2^ 
1 



[(e- £o ) 2 + rV4] [( £ - £o - Wo ) 2 + rV4] 



M 2 T L T R 

+ r 2 ) \u Q 



in 



(eV/2 ~ e ) 2 + T 2 /4 \ 
(eV/2-eo-Lo ) 2 +T 2 /Aj 



In 

2 

T 



(-e1//2- £o + ^ ) 2 + r 2 /4 
(-eV/2 - e ) 2 + r 2 /4 



arctan 



-arctan 



-arctan 



eV/2 - £q 

r/2 

-eV/2 - e 

r/2 

-eV/2 - eo + ujq 

r/2 



arctan 



eV/2 -E -u Q 

r/2 



(CI) 



At zero temperature, we can get an analytical expres- 
sion for heat generation in a single-atom structure by 
using the bare Green's functions in Eq. (|B7i IB9p . Wc 
only take into account the imaginary part of the lead self- 
energies and ignore their energy dependence (the wide- 
band limit )24. Finally, we assume that the phonons are 



uj is the phonon energy, e is the electron onsite energy, 
V is the applied bias, and T = Tl + Tr. The heat genera- 
tion is zero when eV < ujq. Equation (|C1|) can reproduce 
most the qualitative features of heat generation in a sin- 
gle atom, except that it does not take into account heat 
conduction in the phonon system. 
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